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The fractal distribution of charged particles is considered. An example of this 
distribution is the charged particles that are distributed over fractal. The fractional 
integrals are used to describe fractal distribution. These integrals are considered as 
approximations of integrals on fractals. Typical turbulent media could be of a fractal 
structure and the corresponding equations should be changed to include the fractal 
features of the media. The magnetohydrodynamics equations for fractal media are 
derived from the fractional generalization of integral Maxwell equations and integral 
hydrodynamics (balance) equations. Possible equilibrium states for these equations 
are considered. 
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I. INTRODUCTION 

The theory of integrals and derivatives of noninteger order goes back to Leibniz, Liouville, 
Riemann, Grunwald, and Letnikov [l|, y] . Fractional analysis has found many applications in 
recent studies in mechanics and physics. The interest in fractional integrals and derivatives 
has been growing continually during the last few years because of numerous applications. 
In a fairly short period of time the list of such applications becomes long, and include 



chaotic d yna mics 



media 
physics 



18 



_^, 3, 4|, material sciences 15 

1 n r 

11, 1121 , quantum mechanics [13 



2l| , electromagnetic theory 



23l . |24| |. non-Hamiltonian mechanics 125 



diffusion, and transport theory 



yy, 




B|, mechanics of fractal and complex 

2l[ |. astrophysics ^], long-range dissipation 
long-range interaction [271, |28|, anomalous 
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The new type of problem has increased rapidly in areas in which the fractal features 
of a process or the medium impose the necessity of using nontraditional tools in "regular" 
smooth physical equations. To use fractional derivatives and fractional integrals for fractal 



distribution, we must use some continuous model 



10 



ll| . We propose to describe the fractal 



medium by a fractional continuous model lOJ] , where all characteristics and fields are defined 
everywhere in the volume, but they follow some generalized equations that are derived by 
using fractional integrals. In many problems the real fractal structure of the medium can be 
disregarded and the fractal medium can be replaced by some fractional continuous math- 
ematical model. Smoothing of microscopic characteristics over the physically infinitesimal 



iQi,y 



volume transforms the initial fractal medium into the fractional continuous model 
that uses the fractional integrals. The order of the fractional integral is equal to the frac- 
tal dimension of distribution. The fractional integrals allow us to take into account the 
fractality of the distribution. Fractional integrals can be considered as approximations of 



integrals on fractals [38 



391] . In Ref . 38|, l39|] , authors proved that integrals on fractals can 



be approximated by fractional integrals. In Ref. 251], we proved that fractional integrals can 
be considered as integrals over the space with a fractional dimension up to the numerical 
factor. 

The distribution on the fractal can be described by a fractional continuous medium 



Mill 



model |lOl . I 111 . Il2l . |2l|. In the general case, the fractal medium cannot be considered as a 
continuous medium. There are points and domains without particles. In Refs. [lO|, lul, l2l| , 
we suggest considering the fractal distributions as special (fractional) continuous media. We 
use the procedure of replacement of the distribution with fractal mass dimension by some 
continuous model that uses fractional integrals. This procedure is a fractional generalization 
of the Christensen approach 4l| . The suggested procedure leads to the fractional integration 
to describe the fractal medium. In this paper, we consider the magnetohydrodynamics 
equations for the fractal distribution of charged particles. Note that typical turbulent media 
could be of a fractal structure and the corresponding equations should be changed to include 
the fractal features of the media. 

In Sec. II, a brief review of the HausdorfF measure, Hausdorff dimension and integration 
on fractals is suggested to fix notation and provide a convenient reference. The connection 
integration on fractals and fractional integration is discussed. In Sec. Ill, a brief review 
of electrodynamics of the fractal distribution of charged particles is given. The densities 



of electric charge and current for the fractal distribution are described. A fractional gen- 
eralization of the integral Maxwell equation is suggested. In Sec. IV, a brief review of the 
hydrodynamics of fractal media is considered to fix notation and provide a convenient refer- 
ence. In Sec. V, the magnetohydrodynamics equations for the fractal distribution of charged 
particles are derived. The stationary states for these equations are considered. Finally, a 
short conclusion is given in Sec. VI. 

II. INTEGRATION ON FRACTAL AND FRACTIONAL INTEGRATION 

Fractals are measurable metric sets with a noninteger Hausdorff dimension. The main 
property of the fractal is noninteger Hausdorff dimension. Let us consider a brief review of 

nn 

the Hausdorff measure and the Hausdorff dimension [35|, |40| to fix the notation and provide 
a convenient reference. 

A. Hausdorff measure and Hausdorff dimension 

Consider a measurable metric set (W, fin)- The elements of W are denoted hj x,y,z, ..., 
and represented by n-tuples of real numbers, x = (xi, X2, ...,a;„), such that W is embedded 
in M". The set W is restricted by the following conditions: (1) W is closed; (2) W is 
unbounded; (3) W is regular (homogeneous, uniform) with its points randomly distributed. 

The metric d{x, y) as a function of two points x and y & W can be defined by 

n 

d{x,y) = ^\yi-Xi\, (1) 

or 

/ n \ 1/2 

d{x,y) = \x-y\= i^{yi-Xif\ . (2) 

The diameter of a subset E G W cM."' is 

d{E) = diam{E) = sup{d{x,y) : x,y E E}, 

Let us consider a set {Ei} of non-empty subsets Ei such that dim,{Ei) < e, Wi, and 
W C Ui^i ^i- Then, we define 

aE,,D) = uj{D)[dtam{E,)]'' = u;{D)[d{E,)]''. (3) 



The factor uj{D) depends on the geometry of Ei, used for covering W. If {Ei} is the set of 
all (closed or open) balls in W, then 

rrD/22-D 



'^^^) TiD/2 + iy 
The Hausdorff dimension D of a subset E C W is defined 

D = dimH{E) = supjrf G R : 1^h{,E, d) = cxd}, 



32, 


d, 


34 



35|by 



or 



D = diniH{E) = mf{d e R : finiE, d) = 0}. 



From ([5]) and ([6]), we obtain 

1) finiE, d) = ior d> D = dimniE); 

2) fJ^niE, d) = oo for d < D = diinH{E). 



The Hausdorff measure jjlh of a subset E <zW is defined [32|, l33|, [34, l35| by 

oo 

^ih{E,D) = Mm in{{Y^^{E.,,D): E (l[JE,, d{E,) < e Vz}, 



i=l 



or 



/iH(i^, /^) =^P)lim inf V[rf(E,)]^. 



If E C VT and A > 0, then 



where XE = {Ax, x G -E}. 



/iH(A^,/^) = A>H(^,/^ 



(4) 



(5) 



(6) 



(7) 



(8) 



B. Function and integrals on fractal 



Let us consider the functions on W: 

oo 



X) 



i=l 



where xe is the characteristic function of E: Xe{x) = 1 ii x E E, and Xe{x) 
For continuous function /(x): 

lim f{x) = f{y) 



(9) 
Oiix ^ E. 

(10) 



whenever 

\imd{x,y)=0. (11) 

The Lebesgue-Stiehjes integral for ([9]) is defined by 

„ oo 

/ fd^i = y2f3^ME^). (12) 



Therefore 



[ f{x)dfiH{x)= hm \2f{x,)aE„D) 



fW -^-^ - E^ 

It is always possible to divide M" into parallelepipeds: 

Ei^...i„ = {{xi, ...,Xn) eW:Xj = {ij - l)Axj + aj, < aj < Axj, j = 1, ...,n}. (14) 

Then 

dfinix) = lim J{Ei^,„i^,D) = 

n n 

= lim Y\(Axif/" = Y\d''/^''xi. (15) 

The range of integration W may also be parametrized by polar coordinates with r = d{x, 0) 
and angle Q. Then Er^n can be thought of as spherically symmetric covering around a center 
at the origin. In the limit, the function C,{Er^n, D) gives 

dfiHir,n)= lim ^{Er^n, D) = dil^-^r^-^dr. (16) 

diEr,n)^0 

Let us consider f{x) that is symmetric with respect to some center xq G W, i.e., f{x) = 
const for all x, such that d{x, Xq) = r for arbitrary values of r. Then the transformation 

W^W: x^x=x-Xo (17) 



can be performed to shift the center of symmetry. Since W is not a linear space, ( jTTjl 
need not be a map of W onto itself; (ITTj) is measure preserving. Then the integral over a 
D- dimensional metric space is 

/• /"OO 

/ /d/i^ = XiD) / /(r)r^-icir, (18) 

Jw Jo 



where 

^(^' - mry „ ('"' 

This integral is known in the theory of the fractional calculus [1]. The right Riemann- 
Liouville fractional integral is 



(I'^m^) = fp) / (^ - ^r-'mdx. (20) 



For z = 0, Eq. ([20]) gives 



CO 



(/^/)(0) = ^;^y^ x^-'f{x)dx, (21) 



and Eq. ( TTSl) is reproduced by 

Equation (1221) connects the integral on fractal with an integral of fractional order. This result 
permits us to apply different tools of the fractional calculus [Ij for the fractal medium. As 
a result, the fractional integral can be considered as an integral on the fractal up to the 
numerical factor T{D/2)/[2n^/'^T{D)]. 



Note t 
mension 



lat the interpretation of fractional integration is connected with a fractional di- 



251]. This interpretation follows from the well-known formulas for dimensional 



regularization 37|]. The fractional integral can be considered as an integral in the frac- 



tional dimension space up to the numerical factor T {D /2) /[2n^^'^T (D)]. In Ref. 38| it 
was proved that the fractal space-time approach is technically identical to the dimensional 
regularization. 

C. Properties of integrals 

The integral defined in flT3|) satisfies the following properties: 

(1) Linearity: 

/ (a/i + 6/2) dfiH = a fi dfiH + & / /2 dfiH, (23) 

Jw Jw Jw 

where /i and /2 are arbitrary functions; a and b are arbitrary constants. 

(2) Translational invariance: 

/ f{x + Xo)dnH{x) = / f{x)dfiH{x) (24) 

Jw Jw 



since d^H^x — Xq) = dfiH^x) as a consequence of homogeneity (uniformity). 
(3) Scaling property: 

/ f{Xx)dfiHix) = A-^ / fix)dfiH{x) (25) 

Jw Jw 

since dfiH^x/X) = X^^dfi^ix). 

By evaluating the integral of the function f{x) = exp(— ax^ + bx), it has been shown 

37( 1 that conditions (I23|) - (l25l) define the integral up to normalization: 

exp(-ax^ + bx)dnH{x) = vr-^/^a"^/^ exp(&V4a)- (26) 



36 



/ 

Jw 



'w 

Note that, for 6 = 0, Eq. ( l26l) is identical to result from ( !22l) . which can be obtained directly 
without conditions fl23l)-fl25l). 



D. Multi- variable integration on fractal 

The integral in ( IT8|) is defined for a single variable. It is only useful for integrating 
spherically symmetric functions. We consider multiple variables by using the product spaces 
and product measures. 

Let us consider a collection of ra = 3 measurable sets (W4,/ifc,-D) with k = 1,2,3, and 
form a Cartesian product of the sets Wk producing W = Wi x W2 x PV3. The definition 
of product measures and the application of the Fubinis theorem provides a measure for the 
product set W = Wi x W2 y~ W^ as 

(/ii X /X2 X /i3)(w^) = ^il{Wl)^x2{W2)^iz{w^). (27) 

Then integration over a function / on P^ is 

f{xi,X2,X3)d{fii X /i2 X /is) = 
W 

= / / / f{xi,X2,X3)dlJ,i{xi)dl^2{x2)dfi3{x3). (28) 

JWi Jw2 Jw-i 
In this form, the single- variable measure from (ITSll may be used for each coordinate Xk, 

which has an associated dimension at'- 

diik{xk) = \{ak)\xk\'''-^dxk, A; = l,2,3. (29) 



Then the total dimension oi W = Wi x W2 x W3 is 

D = ai + a2 + 03. (30) 

Let us reproduce the resuh for the single- variable integration (TTSl) . from Wi x W2 x W3. We 
take a spherically symmetric function f{xi,X2,X3) = /(r), where r^ = {xiY + {X2Y + {x^Y 
and perform the integration in spherical coordinates {r,(j),6). Equation (l28l) becomes 

d^l{xi)dfi2{x2)dfi3{x3)f{xi,X2, X3) = 

w 



= ^4(0;) / dxi / dx2 / t/xslxil"^ |a;2|"^ ksl"^ f{xi,X2,X3) = 
Jwi Jw2 Jw'i 

= A{a) f dr f d(P /"d^ J3r°^+"2+°^-^(cos0)"^-^(sin0)"2+"^-2(sin^)°^-V(r), (31) 

where J3 = r^sin0 is the Jacobian of the coordinate change, and A{a) = A(ai)A(a;2)A(a3). 

Since the function is only dependent on the radial variable and not the angular variables, 

we can use 

r.n>-.co.-.d.^^(f>r'':{^), (32) 

Jo 2r({ii+i^)/2) ^ ' 

where /x > 0, i/ > 0. From ( l30l) . we obtain 

di^i{xi)di^2{x2)dfi3{x3)f{r) = X{D) / f{r)r^'^dr, (33) 

w J 

where 

MD) ^ ^(^. (34) 



This equation describes the D-dimensional integration 37| of a spherically symmetric func- 
tion, and reproduces the result (fT8|) . 

E. Density function and mass on fractal 

Let us consider the mass that is distributed on the measurable metric set W with the 
fractional Hausdorff dimension D. Suppose that the density of mass distribution is described 
by the function p(r) that is defined by ([9]). In this case, the mass can be derived by 

MoiW) = [ pir)dVD, (35) 

Jw 

where 

dVo = dfii{xi)dfi2{x2)dfi3{x3) = C3{D,r)dxdydz, (36) 



C3(D, r) = A(ai)A(a2)A(a3)x"^-'y"^-'^"-^-\ (37) 

dimHiW) = D = ai + a2 + a^. (38) 

As a result, we have 

Mn{W) = [ p{r)dVD, dVo = C3{D,r)dV,, (39) 

Jw 

where dV^ = dxdydz for Cartesian coordinates, and 

r(ai)r(a2)r(a3) 

As a result, we get the Riemann-Liouville fractional integral [1|] up to a numerical factor 

87r^/2. 



F. Mass of fractal distribution 

The cornerstone of fractals is the noninteger dimension. The fractal dimension can be 
best calculated by box counting method, which means drawing a box of size R and counting 
the mass inside. This mass fractal dimension can be easy measured for fractal media. The 
properties of the fractal medium like mass obey a power law relation M ~ i?^, where M 
is the mass of the fractal medium, i? is a box size (or a sphere radius), and D is a mass 
fractal dimension. The power law relation M ~ R^ can be naturally derived by using 
the fractional integral 10|. The mass fractal dimension is connected 10[ with the order of 
fractional integrals. 

Consider the region W in three-dimensional Euclidean space M'^. The volume of the 
region W is denoted by Vd{W). The mass of the region W in the fractal medium is denoted 
by M£,{W). The fractality of the medium means than the mass of this medium in any region 
W of Euclidean space M? increases more slowly than the volume of this region. For the ball 
region of the fractal medium, this property can be described by the power law M ~ i?^, 
where R is the radius of the ball W . 

The fractal medium is called a homogeneous one if the power law M ~ R^ does not 
depend on the translation of the region. The homogeneity property of the medium can 
be formulated in the form: For all regions W and W of the homogeneous fractal medium 
with the equal volumes VoiW) = Vd{W'), the masses of these regions are equal MuiW) = 
Md{W'). Note that the wide class of the fractal media satisfies the homogeneous property. 
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In Refs. [10(, the continuous medium model for the fractal media was suggested. The 
fractality and homogeneity properties can be realized in the following forms: (1) Homogene- 
ity: The local density of homogeneous fractal media is a translation invariant value that 
has the form p(r) = po = const. (2) Fractality: The mass of the ball region W oi a fractal 
medium obeys a power law relation, M ~ R^, where < D < 3, and R is the radius of 
the ball. These requirements can be realized by the fractional generalization (139|) of the 
equation 

M^iW) = [ p{r)dV,. (41) 

Jw 

The form of function c^lD, r) is defined by the properties of the fractal medium. Note that 
the final equations that relate the physical variables have a form that is independent of a 
numerical factor in the function C3{D,r). However, the dependence of r is important to 
these equations. 

Equation ( l39l) describes the mass that is distributed in the volume and has the mass 
fractal dimension D by fractional integrals. There are many different definitions of fractional 

n 

integrals [l]. The fractional integrals can be used to describe fields that are defined on the 
set W with fractional Hausdorff dimension dimH{W) = D. 
For the Riemann-Liouville fractional integral, 

l/rlai — l|y|«2 — 1 Inlets — 1 

'''^■■->- r(„.)rL)rM ■ (^^' 

where x, y, z are Cartesian's coordinates, and D = ai + a2 + as, < D < 3. 

Note that for D = 2, we have the fractal mass distribution in the volume. In general, 
this case is not equivalent to the distribution on the two-dimensional surface. 

For p(r) = p(|r|) , we can use the Riesz definition of the fractional integrals 1|, and 

csiD,v) = \{D)\vf~^ (43) 

where 

A(«)=7.-'(B)=54;t^. (44) 



Note that 



l\m_^,\D) = {An'/')-\ (45) 



Therefore, we suggest using 

T{D/2) 



X{D) = {4n'/'H\D) = ^^-^^^. (46) 
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The factor (H6|) allows us to derive the usual integral in the limit D -^ (3 — 0). Note that 
the final equations that relate mass, moment of inertia, and radius are independent of the 
numerical factor \{D). 

For the homogeneous medium (p(r) = po = const) and the ball region W = {r : |r| < R}, 



M.(W)^,„i^^£|r|-,v.3. 



Using the spherical coordinates, we get 

r2^-^r(3/2) 



PoR""- 



2^-^TiT{3/2) rj 



DT{D/2) 

As a result, we have M{W) ~ R^, i.e., we derive the equation M ~ R^ up to the numerical 
factor. Therefore the fractal medium with noninteger mass dimension D can be described 
by fractional integral of order D. 

III. ELECTRODYNAMICS OF FRACTAL DISTRIBUTION OF CHARGED 

PARTICLES 

In this section, a brief review of electrodynamics of fractal distribution of charged particles 



2l| is considered to fix notation and provide a convenient reference. 



A. Electric charge for fractal distribution 

Let us consider charged particles that are distributed with a constant density over a 
fractal with Hausdorff dimension D. In this case, the electric charge Q satisfies the scaling 
law Q{R) ~ R^ , whereas for a regular n-dimensional Euclidean object we have Q{R) ~ -R"". 

The total charge of region W is 

Q3m= I p{r,t)dV3, (47) 

Jw 

where p(r, t) is a charge density in the region W. The fractional generalization of ( 1471) is 

Qd{W)= [ p{r,t)dVD, 
Jw 
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where D is a fractal dimension of the distribution, and 

dVD = c,{D,r)dV3. (48) 

The functions c^lD, r) is defined by the properties of the distribution. 

If we consider the ball region W = {r : |r| < R}, and spherically symmetric distribution 
of charged particles (p(r, t) = p{r,t)), then 

_23-^r(3/2) [^ 



«-(«'-^'^T(W^/^W^°"*- 



For the homogeneous case, p{r,t) = po, and 

The distribution of charged particles is called a homogeneous one if all regions W and W 
with the equal volumes Vd(W) = VjjiW') have the equal total charges on these regions, 
Qd{W)=Qd{W'). 

B. Electric current of fractal distribution 

For charged particles with density p{r,t) flowing with velocity u = u(r, t), the current 
density J(r, t) is 

J(r,t) =p(r,t)u. 

The electric current I{S) is defined as the flux of electric charge. Measuring the field J(r, t) 
passing through a surface S = dW gives 

I{S) = ^jiS)= /(J,rfS2), (49) 

Js 

where dS2 = dS2'n is a differential unit of area pointing perpendicular to the surface S, and 

the vector n = rifce/c is a vector of normal. The fractional generalization of fH^ is 



1(3)= /(J(r,t),rfSrf), 

where 

dS, = C2{d,r)dS2, C2{d,r) = A^|r|^-2^ (50) 

Note that C2(2,r) = 1 for d = 2. The boundary dW has the dimension d. In general, the 
dimension d is not equal to 2 and is not equal to {D — 1). 
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C. Charge conservation for fractal distribution 

The electric charge has a fundamental property established by numerous experiments: 
the velocity of charge change in region W bounded by the surface S = dW is equal to the 
flux of charge through this surface. This is known as the law of charge conservation: 

dQ{W) 
dt 

or, in the form 



-I{S), 



^ / p{v,t)dVo = -i (J(r,t),rfS,). (51) 

"f^ Jw JdW 

In particular, when the surface S = dW is fixed, we can write 

Using the fractional generalization of the Gauss's theorem (see the Appendix), we get 

(J(r,t),rfS,) = 
aw 

= j c^\D,r)-^^(c2{d,r)Jk{r,t)yVD. (53) 

The substitution of Eqs. (!52|) and (!53|) into Eq. (ISTjl gives 

C3(Ar)^^ + ^(c2(rf,r)J,(r,t)) =0. (54) 

As a result, we obtain the law of charge conservation in differential form (154|) . This equation 
can be considered as a continuity equation for fractal distribution of particles 11|. 



D. Electric field and Coulomb's law 

For a point charge Q at position r', the electric field at a point r is defined by 

Q v-v' 

Atteo |r — r I'^ 

where Eq is a fundamental constant called the permittivity of free space. 
For a continuous stationary distribution p(r'). 



14 
For Cartesian's coordinates dV^ = dx'dy'dz'. The fractional generalization of (155|) is 

where dV[) = C3{D,r')dV^. Equation (156!) can be considered as Coulomb's law for a fractal 
stationary distribution of electric charges. 

The electric field passing through a surface S = dW gives the electric flux 

^e{S)= /(E,rfS2), 

where E is the electric field vector, and dS2 is a differential unit of area pointing perpendic- 
ular to the surface S. 

E. Gauss's law for fractal distribution 

Gauss's law tells us that the total flux ^e{S) of the electric field E through a closed 
surface S = dW is proportional to the total electric charge Q{W) inside the surface: 

(^E{dW) = -Q{W). (57) 

For the fractal distribution, Gauss's law (1571) states 

E,dS2) = - [ p{r,t)dVD, (58) 

S ^0 Jw 

where E = E(r, t) is the electric field, and p(r, t) is the charge density, dVo = C3{D,r)dV3, 

and ^o is the permittivity of free space. 

If p(r, t) = p{r), and W = {r : \r\ < R}, then 



Q{W)=47i f p{r)c3{D,rydr, 
Jo 



where 

2- 

rJD r^ = - 

T{D/2) 



C3(Ar) = ^^-^SS(^|r|^-^- (59) 



Then 



_23-^r(3/2) r„,^, n_i 



QiW) = 47r ^^^'^^^ ' J^ p{r)r^-^dr. (60) 



For the sphere S = dW = {r : |r| = R}, 

^E{dW)=47rR^E{R). (61) 



15 



Substituting ([60]) and ([6l]) in ([57D, we get 



^w-S^fD^f""^""*- 



E{R) = ^ 2^'''r(3/2) ^^_, ^ ^o_. 



For homogeneous (p(r) = p) distribution, 

23-D 

F. Magnetic field and Biot-Savart law 

The Biot-Savart law relates magnetic fields to the currents that are their sources. For a 
continuous distribution, the law is 

= ^ / ^^^f^dVi, (62) 

47r J^ |r-r'|^ 

where [ , ] is a vector product, J is the current density, and jUq is the permeability of free 
space. The fractional generalization of Eq. ( l62l) is 

This equation is the Biot-Savart law written for a steady current with a fractal distribution 
of electric charges. The law flB5]) can be used to find the magnetic field produced by any 
fractal distribution of steady currents. 

G. Ampere's law for fractal distribution 

The magnetic field in space around an electric current is proportional to the electric 
current that serves as its source. In the case of a static electric field, the line integral of the 
magnetic field around a closed loop is proportional to the electric current fiowing through 
the loop. Ampere's law is equivalent to the steady state of the integral Maxwell equation in 
free space, and relates the spatially varying magnetic field B(r) to the current density J(r). 

Note that, as mentioned in Ref. yjj, Liouville, who was one of the pioneers in the 
development of fractional calculus, was inspired by the problem of fundamental force law in 
Ampbre's electrodynamics and used fractional differential equations in that problem. 

Ampere's law states that the line integral of the magnetic field B along the closed path 
L around a current given in MKS by 
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where dl is the differential length element. For the distribution of particles on the fractal, 

1(8)= /(J,rfS,), 
Js 

where dS^ = C2{d,r)dS2. For the cylindrically symmetric distribution, 

I{S) = 2n f J{r)c2id,r)rdr, 
Jo 

where C2{d,r) is defined in Eq. ( I50ll . i.e.. 

For the circle L = dW = {r : |r| = R}, we get 



f(B,dl) = 2nRB{R). 



As a result, 

B{R) = J^^^-X [ J{ry''dr. 
^ ' RT{d/2) Jo ^ ' 

For the homogeneous distribution, J{r) = Jq, and 

H. Fractional integral Maxwell equations 



Let us consider the fractional integral Maxwell equations 2l|. The Maxwell equations 
are the set of fundamental equations for electric and magnetic fields. The equations that 
can be expressed in integral form are known as Gauss's law, Faraday's law, the absence of 
magnetic monopoles, and Ampere's law with displacement current. In MKS, these become 

E,dS2) = - / pdVn, 
S ^0 Jw 

E,dli) = -^ /"(B,dS2), 

B,dS2)=0, 

s 

B, d\i) = 1^0 I (J, dSd) + ^0/^0 oT / (E, ^82). 
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Let us consider the fields that are defined on the fractal 40] onl y. The hydrodynamic 
and thermodynamics fields can be defined in the fractal media ll|, [l^ . Suppose that the 
electromagnetic field is defined on the fractal as an approximation of some real case with a 
fractal medium. If the electric and magnetic fields are defined on a fractal and does not exist 
outside of the fractal in Eucledian space M.^, then we must use the fractional generalization 
of the integral Maxwell equations in the form 2l| : 



'E,dS, 



1 
^0 Jw 



pdVn, 



[E,d\,) 



|/;b,.s.) 



:B,t/s. 



0, 



(b (B, d\^) = /io / (J, dSd) + eoyuo^ / (E, dSd). 



(64) 



Note that fractional integrals are considered as an approximation of integrals on fractals 

38i,y. 

Using the fractional generalization of Stokes's and Gauss's theorems (see the Appendix), 
we can rewrite Eqs. ( !M|1 in the form 

.-1,^ ^^.^^ 1 



C3\D,r)div{c2id,r)'E)dVD = — I pdVo, 
W ^0 Jw 

J c^\d,T){curl{ci{^,r)E),dSd) = -^ /"(B,dS,), 

/ c^\D, r)div{c2{d, r)^)dVd = 0, 
Jw 

C2^{d,r){curl{ci{-^,r)B),dSd) = Po {3,dSd) + Sopo-^ / {'E,dSd). 
As a result, we obtain 



o?iu(c2(c?,r)E) = — C3(D,r)p, 

c'uH(^Ci(7,r)Ej = -C2(rf,r)— B, 
divic2{d, r)B ) = 0, 



curl (cii'j, r)B j = fioC2{d, r) J + eofioC2{d, r 



dE 
'~dt' 
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Note that the law of absence of magnetic monopoles for the fractal leads us to the equation 

div{c2{d,r)'B) = 0. It can be rewritten as 

(izfB = —(B, gradc2{d,r)). 

In general {d 7^ 2), the vector grad {c2{d,r)) is not equal to zero and the magnetic field 
satisfies divB 7^ 0. If d = 2, we have div(B) 7^ only for nonsolenoidal field B. Therefore 
the magnetic field on the fractal is similar to the nonsolenoidal field. As a result, the 
magnetic field on the fractal can be considered as a field with some "fractional magnetic 
monopole", g^ ~ (B, VC2). 

IV. HYDRODYNAMICS OF FRACTAL MEDIA. 

A. Euler equations for fractal media 



In Ref. Uj, we derive the fractional generalizations of integral balance equations for 
fractal media. These equations leads to the following differential equations. 

(1) The equation of continuity, 

(2) The equation of balance of density of momentum, 

p[^)^Uk = pfk + Vf'pki. (66) 



(3) The equation of balance of density of energy 

d\ / r^ , T-,\ 9'^k 

dxi 
Here, we mean the sum on the repeated index, k and / from 1 to 3, and use the notations 



p(-)^e^c(Z.,.,i.)p -. (67) 



V^A = a{D, d)R'-''-^ (^R^-^A^ . (68) 

= | + a(B.<i)fl«-%A. (69) 

where 

c{D, d, R) = a{D, ci)i?^+i-^, (70) 
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The equations of balance are a set of five equations, which are not closed. These equations, 
in addition to the hydrodynamic fields p(R, t), m(R, t), e(R, t), include also the tensor of 
viscous stress pki(Jl,t). Let us consider the special cases of (l65l) -( l67|) with 



Pki = -pS, 



kh 



where p = p(R, t) is the pressure. Then the hydrodynamics equations fl65|) - fl67|) are 

These equations are the Euler equations for the fractal medium. 

B. Equilibrium equation for fractal distribution 

The equilibrium state of medium means that 

dA _ 9A _ 

dt ' dxk 

for A = {p,Uk,e}. In this case, Eqs. (!72|) - (!74l) give 

fk = ^V^p. (75) 

Equation ( 1751) gives the fractional generalization of the equilibrium equations. From ( l68l) . 
fQlj and ^, Eq. ^ is 



d{c2id,R)p) _ 

^ - PC3[^, rC)fk. 

For the homogeneous medium p{x) = const, and 

d{c2{d,R)p/po) 



C3iD,R)f, 



k 



dxk 
Ifc3(D,i?)/fc = -9f//9xfc,then 

C2((i, R)p + pqU = const. (76) 

This equation is a fractional generalization of the equilibrium equation. 
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C. Fractional Bernoulli integral 



Let us consider Eq. (1731) . Using Eq. (166|) and 



at/D 2 \dt/ D 

we get 

|)^y = W.-^«.Vfp. (77) 

If 

a^/9t = 0, dp/dt = 0, 
then 

Suppose 

/fc = -c(D,rf,i?)9f//9xfc. (79) 

If D = 3 and d = 2, then this force is potential. Using Eqs. ( 178|) and ( 1791) . we get Eq. ( |7711 
in the form 

1(^ + c + pm).o. 

where 

^ ^ r c;(c2(^,i?)p) 

^ ' 4 C2(rf,i?)p • 

As a result, we obtain 

3 2 

XI "^ + ^ + ^^^) = ^^^'^- (^°) 

fc=i 

This integral of motion can be considered as a fractional generalization of the Bernoulli 
integral for fractal media. If the forces fk are potential, and D y^ 3, then the fractional 
analog of the Bernoulli integral does not exist. 
For the density 

p = Poc^\d,R)=po^^R''', (81) 



the integral (IHOl) gives 



2 
For Uk = 0, Eq. §2) leads to Eq. ((761). 



^°" + PoU + C2{d, R)p = const. (82) 
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D. Sound Avaves in fractal media 

Let us consider the small perturbations of Eqs. (!72|) and (173|l : 

P = Po + p', P = Po+ P', Uk = u'k, (83) 

where p' <ti po, and p' <^ po, and po and po describe the steady state: 

^ = ^ = — = — = 
dt ' dxk ' dt ' (9xfc 

Supposing fk = 0, and substituting (183!) into Eqs. (1721) and (1731) . we get 

^ = -PoV^l, (84) 

^ = --Vfp'. (85) 

dt Po 

To derive the independent equations for p', we consider the partial derivative of Eq. ( l84l) 
with respect to time: 

The substitution of (l85l) into (l86l) obtains 



2^/ 



aV' 



For adiabatic processes p = p{p, s), the first order of perturbation is 



V^V^p'. (87) 



p' = v^p', 



where 



As a result, we obtain 



dp'' 
dp) s 



ay 



dt^ 

These equations describe the waves in the fractal medium. 



^'V^V^p' = 0, (88) 



- ^'V^V^p' = 0. (89) 
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V. MAGNETOHYDRODYNAMICS 



A. Magnetohydrodynamics (MHD) equations 



The hydro dynamic and Maxwell equations for a fractal medium UJ, |2ll] are the following. 

(1) The equation of continuity, 

(2) The equation of balance of density of momentum, 

p(^)^«;^ = p/fc-V^p. (91) 

(3) Faraday's law, 

cmH(ci (7, r)E) = -C2(rf,r)-B. (92) 

(4) The absence of magnetic monopoles 

div ( C2{d, r)B) =0. (93) 

(5) Ampere's law, 

CMr/fci(7,r)BJ = yUoC2((i, r)J, (94) 

where the displacement current is neglected. 
Using the Lorenz force density, 

pf=[J,B], (95) 

we get (|9T1) in the form 

We assume a linear relationship between the J and E*: 

J(r,t) = o-E,(r,t), (97) 

where a is the electric conductivity, and E* is an electric field in the moved coordinate 
system. For |u| ^ c, 

E = E, --[u,B]. (98) 

From (EZD, and (El), we get 

E, = a-i J = T^^curl (ci (7, r)B) . (99) 

cr/ioC2(d,r) V / 
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Substitution of ([98]) into ([92]) gives 

CMr/(^Ci(7,r)E*-Ci(7,r)-[u, B] 



dt 
Substituting ([99]) into (ITOO]) . we have 

ci(7,r) 



C2(rf,r)|-B. (100) 



cur/ ( f- — rcur/ ( ci (7, r)B 

VawnCo a, r V 



.cr/ioC2(rf,r) 

-Ci(7,r)i[u,B]) =-C2Kr)^B. (101) 

Then, 

C2(rf,r)|:B = -curl ( ^lhllL^curlU{^,v)B)) + 

+ CMr/(ci(7,r)-[u,B] j . (102) 

As a result, we obtain magnetohydrodynamics (MHD) equations for a fractal distribution 
of charged particles: 

(1) The equation of continuity, 

(Do^-"^"- (™' 

(2) The equation of balance of density of momentum, 

(3) The absence of magnetic monopoles, 

div(c2id,v)B) =0. (105) 

(4) The Ampere law, 

c-ur/ fci (7, r)Bj = /ioC2((i, r) J. (106) 

(5) The diffusion equation for the magnetic field, 

C2(d,r)|-B = -curl ( ^l^^curlL{^,T)B)] + 

+ CMrn ci (7, r)-[u,B]j . (107) 

We have 11 equations for 11 variables p, p, J, u, and B. 
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B. Equilibrium from MHD equations 

Let us consider the stationary (equilibrium) states for MHD equations. The total time 
derivatives in Eqs. (11041) and (I106p are equal to zero, and 



yUoC2(ci 

The substitution of (llOOp into (IllOp obtains 



V^P=[J,B], 
\—-curl{c,{^,r)B). 



D„ 



v> 



Suppose B = {0,0,-82}. Then 



IJ,oC2{d,r) 



CMr/(ci(7, r)B), B 



cur/(ci(7, r)B j = 

= e^(9y(ci(7, T^)Bz) - ey(9^(ci(7, r)^^), 

and 

cwr/ (ci (7, r)Bj, B = 

As a result, Eq. (llOSp gives 

1 



Vfp 



fioC2{d,r 

V^p = ^ 

^ HoC2{d,T 



■BA{ci{i,t)B,), 
-B,dy{ci{-f,v)B,). 



From the definition of V'^, we have 



^C2{d,r)p=- ''^^/\ i?.^(ci(7,r)i?2), 
ox fioC2[d,r) ox 

^C2{d, r)p = - ''^f;'\ i^.^(ci(7, r)Bz). 
dy f^oC2{d,r) dy 



Using AdB = d{AB) - BdA, we get 



— (c id r) + ^3(^,r)ci(7,r) ^2 
dxy"^ ' /ioC2((i, r) "^ 

= ci(7,r)52-:— I r^^B, 



dx \HoC2{d,T) 



(108) 
(109) 

(110) 



111) 



(112) 



(113) 



(114) 
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C2{d,r)p+ ' B 

oy V yUoC2(rf,r) 

-i(7,r)B.#f^%^flJ, (115) 



9 f ^ I J ^\^ , C3(/^,r)ci(7,r) ^2 



f^y VAioC2(c?, r) 
dy 



^{C2{d,r)p) = (}. (116) 



As a result, we obtain 



. , X , C3(£',r)ci(7,r) „ , .„»^ 

C2{d, y)p H — — - — B^ = const. (117) 

HoC2{d,v) 



This equilibrium equation exists only if 

B. ~ ^^^. (118) 

It is easy to see that we do not have the usual invariants for the fractal distribution of 
charged particles. Therefore equilibrium on the fractal exists for the magnetic field that 
satisfies the power law relation 

B, ~ R'^-^+\ (119) 



For the distribution with an integer Hausdorff dimension, we have the usual relation |42 |. 
The typical turbulent media could be of fractal structure, and the corresponding equations 
should be changed to include the fractal features of the media. Therefore, the equilibrium 
of the fractal turbulent medium exists for the magnetic field with the power law relation 

dnni). 

VI. CONCLUSION 

Typical turbulent media could be of a fractal structure, and the corresponding equations 
should be changed to include the fractal features of the media. Magnetohydrodynamics equa- 
tions for the fractal distribution of charged particles are suggested. The fractional integrals 
are used to describe fractal distribution. These integrals are considered as approximations 
of integrals on fractals. Using the fractional generalization of the integral Maxwell equation 
and the integral balance equations, we derive the magnetohydrodynamics equations. Equi- 
librium states for these equations are discussed. The equilibrium for fractal turbulent media 
can exists if the magnetic field satisfies the power law relation. 
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VII. APPENDIX: FRACTIONAL GAUSS'S THEOREM 

Let us derive the fractional generalization of Gauss's theorem, 

/ {J{r,t),dS2)= f div{Jir,t))dVs, (120) 

JdW Jw 

where the vector J(r, t) = J^e^ is a field, and div{J) = di/dv = dJk/dxk- Here, we mean 
the sum on the repeated index k from 1 to 3. Using 

22-d 



we get 



dSd = C2{d,r)dS2, C2{d,r) = :^——\r\'^ ^, 



[ (J(r,t),rfSrf)= / C2id,r){J{r,t),dS2). 
Jaw Jaw 



law Jaw 

Note that 03(2, r) = 1 for d = 2. Using flT^ . we get 



C2((i, r)(J(r, t),rfS2) = / div{c2{d,r)3{r,t))dVs. 

Jw 

23-^r(3/2),^,^^3 



'aw 
The relation 

dVn = CsiD,r)dV3, c^{D,v 



V{D/2) ' ' 

in the form dV^, = c^^{D, r)dVD allows us to derive the fractional generalization of Gauss's 

theorem: 

3{r,t),dSd) = 
aw 



c^^{D,r)divic2{d,r)J{r,t)] dVo- 



'w 
Analogously, we can get the fractional generalization of Stokes's theorem in the form 



where 



/(E,rfg = /c2-'(d,r)(cMH(ci(7,r)E),dSd), 
Jl Js 



2'-T(l/2) , ,_, 
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